home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / src / exampleCode / MP / kernel_switch / mpy_mips1.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.5 KB  |  66 lines
  1. c 
  2. c Copyright (C) 1993, Silicon Graphics, Inc.
  3. c All Rights Reserved.
  4. c
  5.     subroutine mpy_mips1(a,lda,b,ldb,c,ldc,m,n,p,alpha)
  6. c
  7.     implicit double precision(a-h,o-z)
  8.     real*8 alpha,beta
  9.     dimension a(lda,p),b(ldb,n),c(ldc,n)
  10.     integer p
  11. c
  12.     do  100 j=1,(n-1),2
  13.         do 50 i=1,(m-3),4
  14.             t00 = c(i,j)
  15.             t10 = c(i+1,j)
  16.             t20 = c(i+2,j)
  17.             t30 = c(i+3,j)
  18.             t01 = c(i,j+1)
  19.             t11 = c(i+1,j+1)
  20.             t21 = c(i+2,j+1)
  21.             t31 = c(i+3,j+1)
  22.             do 30 k=1,p
  23.             t00    = t00    + a(i,k)*b(k,j)
  24.             t01    = t01    + a(i,k)*b(k,j+1)
  25.             t10    = t10    + a(i+1,k)*b(k,j)
  26.             t11    = t11    + a(i+1,k)*b(k,j+1)
  27.             t20    = t20    + a(i+2,k)*b(k,j)
  28.             t21    = t21    + a(i+2,k)*b(k,j+1)
  29.             t30    = t30    + a(i+3,k)*b(k,j)
  30.             t31    = t31    + a(i+3,k)*b(k,j+1)
  31.  30    continue
  32.                         c(i,j) = t00
  33.                         c(i,j+1) = t01
  34.                         c(i+1,j) = t10
  35.                         c(i+1,j+1) = t11
  36.                         c(i+2,j) = t20
  37.                         c(i+2,j+1) = t21
  38.                         c(i+3,j) = t30
  39.                         c(i+3,j+1) = t31
  40.  50    continue
  41.         do 80 i=i,m
  42.         t00 = c(i,j)
  43.         t01 = c(i,j+1)
  44.         do k =1,p
  45.             t00    = t00    + a(i,k)*b(k,j)
  46.             t01    = t01    + a(i,k)*b(k,j+1)
  47.         enddo
  48.         c(i,j) = t00
  49.         c(i,j) = t01
  50.  80     continue
  51.     
  52.  100     continue
  53.  
  54.         do j=j,n
  55.                 do k=1,p
  56.                         t00 = b(k,j)
  57.                         do i=1,m
  58.                                 c(i,j) = c(i,j) + a(i,k)*t00
  59.                         enddo
  60.                 enddo
  61.         enddo
  62.  
  63.  
  64.     return
  65.     end
  66.